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ABSTRACT 

Neutron stars with very strong surface magnetic fields have been suggested as the site for the origin of observed 
soft gamma repeaters (SGRs). In this paper we investigate the influence of such strong magnetic fields on the 
properties and internal structure of these strongly magnetized neutron stars (magnetars). We study properties of 
a degenerate equilibrium ideal neutron-proton-electron (npe) gas with and without the effects of the anomalous 
nucleon magnetic moments in a strong magnetic field. The presence of a sufficiently strong magnetic field changes 
the ratio of protons to neutrons as well as the neutron drip density. We also study the appearance of muons as well 
as pion condensation in strong magnetic fields. We discuss the possibility that boson condensation in the interior 
of magnetars might be a source of SGRs. 

Subject headings: stars: interiors — stars: magnetic fields — stars: neutron 
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1. INTRODUCTION 

Among the more than two thousand observed cosmologi- 
cal gamma-ray bursts (GRBs), four recurrent sources, so-called 
soft gamma repeaters (SGRs ), have been i dentified and a fifth 
has probably been observed ( Hurley 2000 ). They are believed 
to be a new class of 7-ray transients separate from the source 
of classi cal GRBs. Observations of SGR 0526-66 (M azets et 
al. 1979) SGR 1806-20 (Murakami et al. 1994]) , and SGR 
1900+14 (|Kouveliotou et al. 1998[) with RXTE, ASCA, and 



BeppoSAX have confirmed the fact that these SGRs are newly 
born neutron stars that have very large surface magnetic fields 
(up to 10 15 G) based upon measurements of the spin-down 
timescal e. Recently, SGR 16 27-41 has also been discovered by 
BATSE (IWoods et al. 1999b. It is estimated that its magnetic 



field could be B > 5 x 10 G. The most recent source is SGR 



1801-23 ( pine et al. 2000| ) observed by Ulysses, BATSE, and 
KONUS-Wind. Such st ars have been named magnetars (D un- 
can & Thompson 1992, [Thompson & Duncan 1995 ). [Note 
however that recently Harding et al. (1999) and Marsden et 
al. (1999) have suggested that if relativistic wind outflow con- 
tinuously dominates the spin-down of SGR 1806-20 and SGR 
1900+14, then the surface dipole field may be too low to be 
consistent with a magnetar model.] 

Magnetars have also been suggested as the site for anoma- 
lous X-ray pulsars (AXPs) (v an Paradijs, Taam, & van den 
Heuvel 1995) such as IE 1841-045 (Kes 73) (Gotthelf, Vasisht, 
& Dotani 1 999), RX 107 20.4-3 125 ([ Haberl et al. 1997] ), and IE 
2259+586 (Rho & Petre 1997). However, as another scenario, 



Chatterjee et al. (1999) proposed that these AXPs should have 
magnetic fields similar to neutron stars inferred for ordinary ra- 
dio pulsars and binary X-ray pulsars in an accretion-powered 
X-ray emission model. 

Whether or not magnetars are the source of SGRs or AXPs, 
as relics of stellar interiors, the study of the magnetic fields 
in and around degenerate stars should give important informa- 
tion on the role such fields play in star formation and stellar 
evolution. Indeed, the origin and evolution of stellar magnetic 
fields remains obscure. As early as Ginzburg (1964) and Wolt- 
jer (1964) it was proposed that the magnetic flux ($5 ~ BR 2 ) 
of a star is conserved during its evolution and subsequent col- 



lapse to form a remnant white dwarf or neutron star. A main 
sequence star with radius R ~ 10 11 cm and surface magnetic 
field B 10- 10 4 G [magnetic A-type stars h ave typical surface 
fields < 10 4 G ( [Shapiro & Teukolsky 1983j )] would thus col- 
lapse to form a white dwarf with R ~ 10 9 cm and B ~ 10 5 - 10 s 
G, or a neutron star with R ~ 10 6 c m and B ~ 10" -10 14 G. 
Indeed, shortly after their discovery ( Hewish et al. 1968[ ) pul- 
sars were identified as rotating neutron stars ( Sold 1968| ) with 
magnetic fields B ~ 10 11 - 10 13 G consistent with magnetic field 
amplification by flux conservation. 

Recently, Thompson and Duncan (1993) have invoked a con- 
vective dynamo mechanism to suggest that the magnetic dipole 
field of young neutron stars could realistically reach values of 
the order of 10 14 - 10 15 G, i.e., 10 2 - 10 3 times stronger than 
ordinary pulsars. Moreover, the internal magnetic field of a star 
may not necessar i ly be reflected in its surface magn etic field 
(|Ruderman 1980|, palloway, Proctor, & Weiss 1977b. There- 



fore, the total strength of internal magnetic fields remains un- 
known. Nevertheless, it is expected that appreciably higher 
magnetic fi elds can exist in the interiors of neutron stars ( Rud- 
erman 1980). 

Ultimately, the allowed internal field strength of a star is 
constrained by the scalar virial theorem (cf. C handrasekhar & 
Fermi 1953, |Shapiro & Teukolsky 1983| ), 

2T + W + 3U + M = 0, 

where T is the total kinetic energy, W is the gravitational po- 
tential energy, n is the internal energy, and M. is the magnetic 
energy. For a star of size R and mass M, this gives a maxi- 
mum interior field strength of B ~ 2 x 10 8 (M /M Q )(R/R & T 2 G. 
For neutron stars with R w 10 km and M w M Q , the maximum 
interior field st rength could thus reach B < 10 18 G (Lerc he & 
Schramm 1977). Numerical studies ( pocquet et al. 1995 ) have 
confirmed that neutron stars with ultrastrong magnetic fields are 
stable up to the order of 10 18 G. They also have found that for 
such values the maximum mass of neutron stars increases by 13 
- 29 % relative to the maximum mass of non-magnetized neu- 
t ron stars. This is simila r to the case of magnetic white dwarfs 
([Suh & Mathews 2000d). 



The strength of the internal magnetic field in a neutron star 
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could, in principle, be constrained by any observable conse- 
quences of a strong magnetic field. For example, ra >id motion 
of neutron stars may be due to anisotropic neutrino emission 



induced by a strong magnetic field (e.g., see lanka 1 998). One 



could also consider the effect of ma g netic fields on t tie therma l 



evolution (Heyl & Hernquist 1997, Baiko & Yako /lev 1999) 



and the maximum mass ( [Vshivtsev & Serebryakovja 1994] ) of 
neutron stars. 
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Recently, Chakrabarty et al. (1995) llave inves 
tigated the gross properties of cold nuclear matter 
magnetic field in the context of a relativistic Harjree 
and have applied their equation of state to obtain the 
masses and radii for magnetic neutron stars. 

Since strong interior magnetic fields modify 
equation of state for degenerate stars, their mass 
tion will be changed relative to that of nonmagnetic 
cently, we have obta ined a revised mass-radiu s relatiqn 
netic white dwarfs (Suh & Mathews 2000a) with 
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of state for electrons in a strong magnetic field. For 
nal magnetic fields of B ~ 4.4 x (10 11 - 10 13 ) G, we 
that both the mass and radius increase disting 
mass-radius relation of some observed magnetic 
may be better fit if strong internal fields are assumec 
In this work, we now extend this study to an inve; 
the effect of magnetic fields on the internal propert 
tron stars as well. If ultrastrong magnetic fields exist 
rior of neutron stars, such fields will primarily affect 
ior of the residual charged particles. Moreover 
from the anomalous magnetic moment (AMM) of tl 
in a strong magnetic field should also be significant 
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Prakash, & Lattimer 2000). In particular, in a stron 
field, complete spin polarization of the neutrons oc> 
the interaction of the neutron magnetic moment wit 
netic field. Therefore, we consider both cases with 
the effects of the AMM. 

Even so, standard internal properties such as 
equation of state, neutron drip, and the threshold den 
particles will be modified by a strong magnetic fiek I 
poses of illustration, we will consider a degenerate 
i nteracting neutron-proton-e lectron (npe) gas in 
( [Shapiro & Teukolsky 1983| ). We find that under 
of charge neutrality and chemical equilibrium, the { 
a sufficiently strong magnetic field changes the ratio 
to neutrons as well as the threshold density for the 
of muons and pion condensation. 

In Sec. 2 we first review the properties of 
scribe the equation of state for a particle gas in an 
netic field. In Sec. 3 we first consider the case 
particle AMM and we derive the proton to neutron 
ideal npe gas for the lowest Landau level analytica 
numerically obtain the neutron appearance density 



concentration, Y p , in magnetic fields. In Sec. 4 we 
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effect of the particle AMM on the equation of state 
abatic index. In Sec. 5 we calculate the muon threshold 
in a magnetic field. We consider pion production 
sation in strongly magnetized neutron-star matter 
Finally, we discuss the possibility that boson condensation in 
the interior of magnetars is a source of SGRs. 



2. EQUATION OF STATE FOR PARTICLES IN A MAGNE' 

2.1. Electrons and muons 

Let us begin by conside ring the properties of an 
an external magnetic field ( f^andau & Lifshitz 1938 
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ICanuto & Chiu 1968[ |Schwinger 1988j ). For 
jh field, the energy states of an electron in a 
■e quantized and its properties are modified ac- 
cyclotron energy becomes comparable to the 
ss energy and the electrons in the excited Lan- 
relativistic. 

estigate these effects for electrons we must first 
equation in an external, static, and homoge- 
field. With a uniform magnetic field B along 
convenient choice of gauge for the vector 
o = and A = (-yB, 0,0), we obtain the elec- 
pns and energy dispersion in a magnetic field 



tie 



(lohnson & Lip >mann 1949, Schwinger 1988). The dispersion 



relation relation 
field is 



for an electron propagating through a magnetic 



F ■ = 



n + 



where n^ ■ 
ber of the Landtu 
tion onto the mi 
c is the speed 
tron momenturr 
the electron 



B e c = m 2 c 3 /eh = 
for 7 e <C 1 and t 
where a ~ 1 / 



137 



1988). 

The main modification 
comes from the 



+ s e z , in which n is the principal quantum num- 
level, si = ±1 /2 is the electron spin projec- 
gnetic field direction, e is the electron charge, 
light, h is Planck's constant, pi is the elec- 
along the z-axis, and m e is the rest mass of 
e, let us use a definition j e = B/B e c where 
4.414 x 10 13 G. In Eq. (0), n e = - (a/4ir)~f e , 
= (a/47r)[ln(2 7e )-(C+|)] 2 + - • •, for le » 1, 
and C = 0.577 is Euler's constant (S chwinger 
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of an electron in a magnetic field 
available density of states for the electrons 
(|Landau & Lifsjhitz 1938). The electron state density in the 



absence of a ma 'netic field, 



is replaced with 



level n"f as 



Ce(e e ) 



'2 2 



c A +m 2 c 4 + 2nceBn e j:] l '^ + m e c^K e , 



,1/2. 



(1) 



d 3 p 
(2^P 



eB 
(2^ 



-dp z 



(2) 



(3) 



f eld. 



This modification affects the thermody- 
) of the electron gas. 
consider a gas of electrons at zero temperature in a 



in a magnetic 
namic propertie ; 
Let us 

magnetic field ([ Jlandford & Hernquist 1982). From Eq. (1) we 
can define the el i 



ctron Fermi energy E e F for an arbitrary Landau 



m 2 c 4 + p e F 2 c 2 + 2hceBn e f ] L / 2 + m e c 2 n 



(4) 



Here p e F denote; the electron Fermi momentum. 

Now we can obtain all of thermodynamic quantities in terms 
of the Fermi en ;rgy, Eq. (||), and the phase space integration, 
Eq. (j^), in a magnetic field. The number density of electrons in 
a magnetic field is then given by 



where e P is defined as e e =E P / m e c 2 and 



= J2J2 J(e e -K e ) 2 -{l+2 le n e f ). 



(5) 



(6) 



n=0 si 



Suh & Mathews 
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The maximum Landau level n e max for a given electron Fermi en- 
ergy e e and magnetic field strength j e is given by 



(e e -K e ) 2 -l I e 
U Z 



>n, 



(7) 



where int [x] means an integer value of the argument x. The 
pressure of an ideal electron gas in a magnetic field is then 



7 e , /m e c . 
/ ' = 4^'"'' (.^) *<"<>• 



(8) 



where 



n=0 si 



(e e -K e )J(e e -K e ) 2 -(l+ 2j e n e f ) 



{e e -K e )+J (e e - Kef - ( 1 + 2j e n e f ) s 



-{\+2 le n e f ) In 
Similarly, the energy density is 

C- e (e e ) = -^2 m eC \ -^~ ) Xe(£e), 



(9) 



(10) 



where 



n=0 si 



(e e - K e ) x j (e e -n e ) 2 -(l+ 2j e n e f ) 



(e e -K e )+J (e e - K e ) 2 - ( 1 + 2j e n e f ) s 



+{\+2 le n e f ) In 
From these, we obtain the energy per electron 

2 Xefe) 



(ID 



E e (e e ) = m e c 



Ce(e e ) 



(12) 



2.2. Protons 

Although the proton mass is much greater than the electron 
mass, magnet ic effects on protons c an be as important as those 
on electrons (Lai & Shapiro 1991). For instance, the proton 
pressure is always much smaller than the electron pressure at 
low density. But, ignoring the influence of the magnetic field 
on protons would lead to the unphysical result of proton pres- 
sure dominance at low density. Therefore, whenever the mag- 
netic field significantly affects the electrons, it also affects the 
protons. 



The energy dispersion relation for protons E ,, for an arbi- 

"/ 

trary Landau level in a magnetic field is: 



E p „ = 

"f 



,fc J + ^[{i + 2 V ;} 1/2 - s f*] 2 



1/2 



, (13) 



where rA = n+ 1 - s 1 ', n is the principal quantum number of the 
Landau level, s P = ± A is the z component of the proton spin, p P 
is the proton momentum along the z-axis, and n P N = (en/m p c)K p , 
with k p = 2.79, is the proton anomalous magnetic moment. 
The proton number density in a magnetic field is : 



1p 



H " = 2n 2 



where e p = Ep /m p c 2 and 



' 2 — in 2 



n=0 £ 



l+2j p n p f -s P K p j p . 



(14) 



(15) 



(16) 



The maximum Landau level n p max for a proton in a magnetic 
field is given by 



(e p + S P z K plp ) 2 -\ I l_ sP 
2 1p \ 2 Z , 



The pressure of a proton gas in a magnetic field is 



(17) 



Pp= A^ m 



pc 2 (^ff * P (e P ), (18) 



where 

w=EE 



n=0 



e pj£ 2 p-™ 2 p -™ 2 p ln i 



< „ + ^/e p -m p ^ 



where 



Note, that as -f e goes to zero, Eqs. (||) and (|^) - ( |l2| ) recover 
exactly the usual non-m agnetic equation of state for electrons 
dSuh & Mathews 2000a| ). 

We can obtain similar quantities for muons simply by re- 
placing the electron quantities by the corresponding muon 
quantities (e.g., replace j e = B/B e c by 7 M = B/B£, where B^ = 
m 2 ^ jeK). The muon AMM k m has nearly the same value as j "*» 

K e . The difference is only k m — n e ~ 0.59 x 1(T 5 (Schwinger Xp(tp) = , E E 

1973) [The exp erimental value for this difference is 0.63 x 10~ 5 
(prandy 199 if)]. 



The energy density is 

c i , IP 2( m P c \ i , x 

£p(£p)= ~^2 m P c {— ) Xp(e P ), 



(19) 



(20) 



n=0 s p - 



e p Je 2 p -m p +ml In 



2.3. Neutrons 



' e P + \h 2 P -™V 



(21) 



A neutral Dirac fermion can interact with an external elec- 
tromagnetic field by means of the Pauli non-minimal coupling. 
Then the energy dispersion relation E F for neutrons in a mag- 
netic field is given by 



E" F = [p 2 c 2 + (Jmy + p 2 ± c 2 + s n z ^B) 2 ] 1/2 (22) 



where s" = ± \ is the neutron spin projection onto the magnetic 
field direction, p» and are the components of the neutron 
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momentum parallel and perpendicular to the magnetic field, and 
H n N = (eH/m n c)K n , with n„ = —1.91, is the neutron anomalous 
magnetic moment. 

In the calculations of thermodynamic quantities, the phase- 
space integration separates into two steps which involve inte- 
gration ov er pii and p \ (e.g., see B roderick, Prakash, & Lat- 
timer 2000). Then the neutron number density in a magnetic 
field is given by 



1 /m„c\ 3 



(23) 



where 



Cn(e«) = ^2 \ 



[^--„ 2 ] 3/2 



+ 2 s z ^n"in 



h n yje 2 n -ml + e 2 r j arcsin ) - | } 



(24) 



In the above Eqs ( ^3| ) and (Q), we have defined e„ = E n F jm n c 2 , 
7„ =B/B n c ,B n c = eh /m 2 n c 3 , and 



m„ = 1 +i"fc„7„. 
The pressure of a neutron gas in a magnetic field is 

3 



where 



1 9 /m n c\ 3 



(25) 

(26) 
(27) 



In equation (27), 

$° = (2 e 2-4-m2 + 9 K „ 2 7 „> + (4-m2-9^ 7 2)^Q, (28) 



2P+4{m2+m„+6(< 7 2-l)}Q 



10 



el + %{K 2 nl 2 n -\) 



where we defined the quantities: 



n, 



(29) 



e„i/ej-mj, Q = ln 



„ I . fnin\ 7T 

ft = e^arcsin(-)-^ 



(30) 



Finally the energy density of a neutron gas in a magnetic field 



S n (e„)= —^rrinC 2 \^-J (31) 



where 



and 



X„ + s z k„7«X„ 



(32) 



Xl = 3 [(2e* - m 2 n )V - m 4 , Q] , = 4 [m„P + m 3 Q + 2efc] 



(33) 



3. INVERSE /3-DECAY AND NEUTRON APPEARANCE IN A STRONG 
MAGNETIC FIELD: WITHOUT THE ANOMALOUS MAGNETIC 
MOMENTS 

Let us first consider the physics of an npe gas in a strong 
magnetic field without the anomalous magnetic moments of 
particles, i.e., Kj = Q(j = e,p,n). This illustrates the domi- 
nant physics at moderate magnetic field strength. For illus- 
tration, consider a homogeneous gas of free neutrons, protons, 
and electrons in /3-equilibrium in a uniform magnetic field. At 
low densities, the most energetically favorable nucleus is 56 Fe 
which is the endpoint of thermonuclear reactions. As the den- 
sity increases above ~ 10 4 gem" 3 , electrons become unbound 
and relativistic. At sufficiently high densities, p > 8 x 10 6 

gem" 3 , protons in nuclei are converted into neutrons via inverse 
/3-decay: 

e~ + p — >n + v. (34) 

Since the neutrinos can escape, energy is transported away 
from the system. Thus, the equation of state in the star will 
be modified mainly due to the inverse /3-decay. The reaction 
([34]) can proceed whenever the electron acquires enough en- 
ergy to exceed the mass difference between protons and neu- 
trons, Q = m„-m p = 1 .293 MeV. The transformation of protons 
into neutrons, reaction (Q), is effective whenever the /3-decay 
reaction; 

n — ► p + e + v (35) 

is slower than the rate of electron capture by protons. Reaction 
( |35| ) is blocked if the density is high enough that all energeti- 
cally available electron energy levels in the Fermi sea are occu- 
pied. Thus, there is a critical density for the onset of reaction 

&■ 

Similar to the field-free case, we can take into account the 
above processes in an intense magnetic field. Assuming that a 
mixture of free neutrons, protons, and electrons are in equilib- 
rium, then reaction (53) implies 



p e + Pp = p n , 



(36) 



where pj = E J F (j = e,p,n,^,7r) is the chemical potential of the 
jth particle. We have set the neutrino chemical potential [i Ve to 
zero. Let us now define 



Pf Mj a \ n 

x; = , e,= — -z, and A,- = 

nijC nijC 1 nijC 



(37) 



where \j is the Compton wavelength of the jth particle. 
From chemical equilibrium, Eq. (Da), we have 



tn e c 2 e e + m p c 2 e p = tn„c 2 e„, 



and charge neutrality gives 



t e — fl p- 



(38) 



(39) 



In order to determine the equilibrium composition and hence 
the equation of state, the above equations should be solved si- 
multaneously. 

Consider now the minimum density at which neutrons first 
appear in a strong magnetic field. This neutron appearance 
density is determined by setting n„ = 0, or e„ = 1 = (1 +x 2 ) 1 ! 2 . 



Suh & Mathews 
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Since the protons at this density are nonrelativistic, i.e., e* « 1, 
we approximately obtain e* ~ 2.53 as the specific electron 
chemical potential at which neutrons first appear according to 
Eq. ([38j) (Hereafter an asterisk is used to denote a threshold 
value for the appearance of new particles.) Therefore, we have 
n e max = for electrons if j e > (e* 2 - l)/2 ~ 2.7. That is, for 
B > 2.7 B e c ~ 1.2 x 10 14 G, electrons reside in the lowest Lan- 
dau level. Substantially n p maK may not be zero for 2^ p <C 1. But, 
in order to compare the two cases of higher-Landau-level and 
lowest-Landau-level occupation for the charged particles, we 
simply take n p maK = for protons. 

Since Q and m e are both much less than m n , from Eqs. ( |3~8| ) 
and ([39]) we then obtain the proton-to-neutron ratio analytically 
when we assume that electrons and protons are in the lowest 
Landau level. Then 



(2m n Q + C 2 „nl /3 ) 2 -4m 2 n m 2 
4C 2 p n 2 n (m 2 n + C 2 n 2 J 2 ) 

where C p = 2n 2 m p X p /^ p , and C 3 = 37r 2 m 3 A 3 . 



1/2 



(40) 
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strong magnetic field through the direct URCA process (B aiko 
& Yakovlev 1999, Leinson & Perez 1998| ). [Note thatBaiko & 
Yakovlev (1999) have considered the direct URCA process at 
the core of a neutron star (p > 10 14 gcm~ 3 ) and not its crust. In 
reality, the direct URCA process can never proceed at such low 
density as 10 12 gem" 3 because most protons are confined within 
the nuclei at p < 10 14 gcm~ 3 in realistic neutron-star matter]. 

However, electrons and protons, actually, are not in the low- 
est Landau level for higher densities. Above a critical density, 
higher Landau levels begin to contribute to the chemical po- 
tential of the electrons and protons and hence particle number 
densities. Ultimately, discrete Landau levels become continu- 
ous and the proton concentration Y p reverts to the nonmagnetic 
limit as the density increases. As a result, inverse (3 decay can 
still be suppressed at high densities in strong magnetic fields. 
Therefore, neutron-star rapid cooling may not be affected by 
the direct URCA process even though it is enhanced in strong 
magnetic fields. However, in order to enhance the cooling by 
the direct URCA process, one can invoke other mechanisms 
such as boson condensation (Tsuruta 1998), nucleon superflu- 
idity (Yakovlev et al 1999), etc., if they exist. 

The proton-to-neutron ratio, Eq. (|40[), gives the number den- 
sity at which neutrons first appear: 



nl(B) = n*<B) = n* e (B) = 



2tt 2 A 3 



Q 1 



1/2 



(41) 



Comparing with the zero-field result 



3n z Xi 



Q 2 -™ 2 e 



-1 3/2 



(42) 



we obtain the relative density at which neutrons appear in a 
strongly mag netized neutron star for B > 2.1 Bj (e.g., see L ai & 
Shapiro 1991): 



p,:qb) _ 3 



2 le 



Q 2 — m 2 



■■0.211 le 



(43) 



Fig. 1 The proton fraction Y p = rip/ng without the AMM as a function 
of the neutron density p„ (in unit gem -3 ) for the given value of log7 f , 
7„ = B/Bc's. The 6 = line is the non-magnetic case. The dashed lines 
occur if charged particles are restricted in the lowest Landau level. 



We can see that the neutron appearance density increases lin- 
early with the magnetic field B. This result is equivalent to one 
direc tly calculated from the general form ( Shapiro & Teukolsky 
1983), 



m p n*{B) 



Figure 1 shows the proton fraction Y p = n p /ftB (where tig = 
«« + ftp) as a function of the neutron density p„ for values of j e 
less than 10 3 . (For higher fields, we should take into account 
the nucleon AMM). In the nonmagnetic case, the conditions of 
charge neutrality and chemical equilibrium in a npe gas lead 
to a threshold for an increase in the proton concentration up to 



a value of Y p ~ 1/9 as p„ exceeds w 10 12 gem 3 . This means 



that the inverse (3 decay is strongly suppressed by Pauli block- 
ing in neutron-rich nuclear matter which consists only of a npe 
gas. However, in the case of a strong magnetic field, if we as- 
sume that electrons and protons are always in the lowest Landau 
level, then an increase in the concentration of protons does not 
occur even as p n exceeds 10 12 gem" 3 . That is, inverse (3 decay 
is not suppressed in magnetic fields. Far from suppressing the 
inverse decay, the magnetic field instead catalyzes the reac- 
tion. This means that rapid neutron-star cooling can occur in a 



e n=0 si 

with e* ~Qlm e . 

4. EFFECTS OF ANOMALOUS MAGNETIC MOMENTS 

Since n n = at neutron appearance, we have e* = 1 +s"k,„^„. 
Then, from chemical equilibrium and charge neutrality, we ob- 
tain the electron Fermi energy at neutron appearance e* when 
the nucleon AMM is included: 



r] 2 -tn 2 „ + m 2 e 



2m e r\ 



(45) 



where 77 = m„(l + s"K„j„) + m p s p k£y p . Note that rj depends on 
both s" and s?. Of course Eq. ( p3| ) goes to e* ~ Q/m e when 
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K Pt „ = 0. In Eq. (^5() rj is similar to the neutron effective mass in 
a magnetic field because 77 becomes m„ when n pn = 0. Hence, 
for = -1/2, it is possible for 77 to have a negative value as 
7 C increases. But this case is unphysical because e* becomes 
negative. Therefore, we take the proton spin s 1 ' = +1 /2 in this 
work. Now we can obtain the neutron appearance density when 
the nucleon AMM is included 



P*(B) _ 3 
PJ(0) 2(e 



*2. 



1)3/2^^ 

7 M=0 J' 



(e*-K e ) 2 -(l+2 7e np. (46) 



Here, n e * ax denotes the maximum electron Landau level at neu- 
tron appearance. The neutron threshold density p*(B) is plotted 
in figure 2 as a function of 7e = B/B e c . We find that p* is signif- 
icantly affected by the nucleon AMM above log 7 ,» w 3. Also, 
above this field strength, the neutron appearance density is split 
according to the orientation of the neutron spin. 




Fig. 2 Neutron appearance density p,* (in gem -3 ) as a function of 7 e . 
The solid line is the case without the AMM. The dashed line corre- 
sponds to s" = -1 /2 when the nucleon AMM's are included. 

Now the equation of state for a npe gas at p > p*(B) can be 
determined in terms of the parameter e e . The charge neutrality 
condition, Eq. (p3[), at a fixed j e , 



£X>£fov# = (47) 



n=0 



n=0 



gives e p for a given e e . We also obtain e„ from the chemical 
equilibrium, Eq. ([32|). Note that both sides of Eq. ( |47| ) have the 
same statistical weight for the excited Landau level when the 
AMM of the proton is ignored. However, the statistical weight 
of protons is changed when the proton AMM is included. 

Finally, we obtain the total baryon density tig = n p +n n and the 
proton concentration Y p = n p /n B . In figure 3 we can see that the 
nucleon AMM above a field strength of j e > 10 3 significantly 
affects the value of Y p . The difference between the AMM and 
non-AMM results above a density of > 10 12 gem" 3 comes 
from the polarization of the proton spin. Besides, the equation 
of state, the mass-energy density p = (£ e + £ p + £„)/c 2 and the 



pressure P = P e +P p +P n , are straightforwardly determined. In 
calculations of the equation of state and the adiabatic index, we 
take s n z =-1/2. 
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Fig. 3 The proton fraction Y p = n p /iiB with the AMM included as a 
function of the neutron density p„ (in unit gcirT 3 ) for the given value 
of log7 e . The solid line is for the case without the nucleon AMM. The 
dashed lines are the cases of s" = -1/2. 



Lai & Shapiro (1991) have shown the equation of state and 
the adiabatic index for a field strength less than log 7e = 2. In 
this work, we calculate those for a field strength less than log7 e 
= 5. Figure 4 shows the equation of state for a strongly mag- 
netized npe gas. At low density (p < p*) and high magnetic 
field (B > 4.4 x 10 14 G), we utilize the fact that electrons are 
nonrelativistic and are in the lowest Landau level as in Lai & 
Shapiro (1991). We plot the adiabatic index Y = dlnP/dlnp as 
a function of p in Figure 5. We find a non-zero Y at the neutron 
appearance as the magnetic field increases. 

The adiabatic index is a crucial factor for understanding the 
global radial stability of a st ar as well as the local sound speed 
(Shapiro & Teukolsky 1983). In figure 5 we know that at neu- 



tron appearance Y has a non-zero value and the minimum value 
of r increases as the magnetic field strength increases. But the 
nucleon AMM reduces the minimum of Y at neutron appear- 
ance. A strange feature in figure 5 is that for a magnetic field 
strength greater than ~ 10 15 G (log7 e > 2), an oscillatory be- 
havior in the adiabatic index begins to appear above a density of 
p > 10 12 gem" 3 . This means that at high densities and fields 
Y is significantly affected by the proton fraction even though it 
is small (Y p < 0.1). Notice that for log7 e < 2, as the density in- 
creases above the neutron appearance density p*, the oscillatory 
behavior vanishes since neutron pressure dominates. 

This oscillatory behavior has the same physical origin as 
the well-known de Haas-van Alphen effect ( Landau & Lifshitz 
1938). It arises as electrons begin to fill the next unoccupied 
Landau level. Increasing the density increases the occupation 
of this level and does not lead to a rapidly increased pressure. It 
would be interesting to understand the physical consequences 
of the oscillatory behavior of the adiabatic index in strongly 
magnetized neutron stars. We speculate that this might drive a 
pulsational instability not unlike classical cepheids. That is, as 
a region develops low Y it may be unstable to collapse. At high 
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density, however, T suddenly stiffens and the region bounces as 
the r is restored to a higher value. To explore this possibility, 
however, it would be necessary to calculate the evolution of the 
interior of strongly magnetized neutron stars in a realistic mag- 
netohy drodynamical model. T his will be the subject of a future 
study (|Suh & Mathews 2000b|). 




Fig. 
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logp 



4 Pressure P (in dyn/cm 2 ) vs. total density p (in gem -3 ) for a 
npe gas in magnetic fields of various strengths, log7<,. The 6 = line 
corresponds to the non-magnetic case. The solid line are for the case 
without the AMM. The dashed lines are for s'l = -1/2 with the AMM 
included. Numbers are labeled as in Fig. 3. 
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FIG. 5 Adiabatic index F vs. p (in gem ) for a npe gas in various 
magnetic fields. Curves are labeled as in Fig. 4. 

5. MUON PRODUCTION IN A STRONG MAGNETIC FIELD 

In order to study the appearance of new particles at high 
density, let us consider muons in an ideal npe gas. Normally, 
muons decay to electrons via, 



But when the Fermi energy of the electrons approaches the 
muon rest mass ~ 105.66 MeV, it becomes energetically 
favorable for electrons at the top level of Fermi sea to decay 
into muons with neutrinos and antineutrinos escaping from the 
star. Hence, above some density, muons and electrons are in 
equilibrium: 



fi «-> e 



(49) 



assuming that the neutrinos leave the star. In this chemical equi- 
librium, charge conservation implies 



Mm = AV 

Equilibrium between n, p, and e means 

ftn — ftp ftei 

and charge neutrality requires 



(50) 



(51) 



(52) 



Now consider the minimum density at which muons are first 
produced in a strong magnetic field. The threshold condition 
for muons to appear is given as n M = 0. In order to satisfy this 
condition, e M must be unity. If 1 , the muon number density 
rift is not zero even though the maximum Landau level «m ax for 
muons in a strong magnetic field is zero. Thus, from Eq. (p0[), 
we simply obtain 

e ; = ^. (53) 



Then e* and e* are given by chemical equilibrium, Eq. (|51|), and 



the charge neutrality condition, Eq. 

n M = 0. Actually, since e* and e* depenc ...^ 

magnetic field, we solve Eqs. ( pl| ) and (|52|) simultaneously to 
obtain 



, for a given 7 C when 
on the strength of the 



and 



/ m 2 „ + m}, — m? 

* IP M e n 

W ml S * Kplp > 



— e + — 

p 

m„ ' m„ 



(54) 



(55) 



Thus, the threshold density for the appearance of muons in a 
magnetic field becomes: 



j=e,p J 



(56) 



(48) 



Figure 6 shows the muon threshold density p* (B) as a function 
of 7 e . We can see that pt(B) without the AMM is not affected 
by magnetic fields less than B ~ 10 18 G. But when we consider 
the AMM, its contribution to p*(B) becomes important above 
B ~ 10 17 G. However, the equation of state for a npe gas with 
or without a magnetic field is nearly unaffected b y the existence 
of muons in neutron stars ( |Canuto & Chiu 1968 ). 

We should als o correct for the appearance of hyperons (B ethe 
& Johnson 1974). Indeed light hyperons with mass less than 
~ 1 .2 GeV are expected to appear at typical neutron star densi- 
ties depending on the model. But the resulting equation of state 
is not very differe nt from that of pure neutron matter ( Shapiro 
& Teukolsky 1983). 
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Fig. 6 Appearance density p* (in gem -3 ) tor muons and pions 
(j = p~ , 7r~) as a function of y e . The solid line are for the cases without 
the AMM. The dashed lines are the cases with the AMM included. 



6. PION PRODUCTION AND CONDENSATION IN A STRONG 
MAGNETIC FIELD 

At very high density (p > p*), neutron-rich nuclear mat- 
ter is no longer the true ground state of neutron-star matter. It 
will quickly decay by weak interactions into chemically equi- 
librated neutron-star matter. Fundamental constituents, besides 
neutrons, may then include a fraction of protons, hyperons, and 
possibly more massive baryons. Moreover, a phase transition 
to quark matter and boson (pion, kaon) condensation are also 
possible. However, a first-order phase transition by thermal 
nucleation to quark matter may not occur in a magnetic field 
(Chakrabarty 1995). Boson condensation in an external mag- 
netic field is also a very subtle problem. 

Recently, a true Bose-Einstein condensation (BEC) has been 
experimently realized in a system of 87 Rb atoms that was con- 
fined by magnetic fields and evaporatively cooled (A nderson et 
al. 1995). An important consequence of the possible appear- 
ance of spin zero bosons is that they can form a Bose-Einstein 
condensation (BEC) at sufficiently low temperatures. An ideal 
condensation consists of a large number of bosons in a state of 
zero kinetic energy. This would have at least two implications. 
One is that the equation of state would be softened, and the 
other is th at the cooling rate from esca ping neutrinos would be 
enhanced (Shapiro & Teukolsky 1983). 

Regarding the BEC in a magnetic field, long ago Schafroth 
(1955) pointed out that for a non-relativistic boson gas, BEC 
does not take place in the presence of a constant magnetic field. 
It also was shown by Toms (1994, 1995) that generally a BEC 
in the presence of a constant magnetic field does not occur in 
any number of spatial dimensions. However, recently Elmfors 
et al. (1995) suggested that condensation may occurs in three 
dimension since the Landau ground state can accommodate a 
large charge density even though it is not exactly a BEC. In 
particular, Rojas (1996) has shown that a BEC actually may oc- 
cur in the presence of a constant homogeneous magnetic field in 
three dimension without requiring the vanishing of the chemical 
potential. There is, however, no critical temperature at which 
condensation begins. 



Actually, since the criterion for condensation (usually taken 
as the equality of the chemical potential to the ground state en- 
ergy) leads to a divergence of the number density, condensa- 
tion cannot occur in a magnetic field. However, if the chemical 
potential depends on both temperature and the magnetic field 
the d ivergence is avoided and condensation may occur (R ojas 
1996). Eventually, in a strong magnetic field and at any non- 
zero temperature, the number density in the ground state be- 
comes finite. Therefore, in this work we will assume that boson 
condensation occurs in a strongly magnetized neutron star and 
that the boson number density has a finite value at sufficiently 
low temperature. 

In order to obtain the energy dispersion relation for a spin 
zero charged boson, we solve the Klein-Gordon equation in 
an external magnetic field. Under the same conditions as for 
charged fermions of Section 2, the dispersion relation for a 
charged boson in a magnetic field is then given by 



E b = [pic 1 + m 2 b c 4 + hceBn b ] 



1/2 



(57) 



where nt, = 2n+ 1, (n = 0, 1, . . .), n is the principal quantum 
number of the Landau level, and b denotes charged bosons 
(b = 71" , K^, ■ ■ ■)■ Note that for charged bosons in a magnetic 
field, their energy state depends on the magnetic field strength 
even though they are in the lowest Landau level (n = 0). Simi- 
larly, the boson state density in the absence of a magnetic field, 



1 



is now replaced with 



n c 

ra=0 



d 3 p 



eB 



:dp z 



(58) 



(59) 



in a magnetic field. 

Neutron stars will provide a unique opportunity to verify 
the hypothesis of boson (pion, kaon) condensation in a strong 
magnetic field. In this work, as an example, we only consider 
charged pion condensation via n — > p + ir~. Since an ideal cold 
npep, gas allows 7r"'s to be produced, let us consider the pion 
appearance threshold in a strongly magnetized neutron star. If 
we neglect the strong interaction between pions and nucleons, 
negatively charged pions are formed through the reaction: 



■ p + ir 



(60) 



in dense neutron matter when the electron chemical potential 
p e exceeds the ir~ rest mass, = 139.6 MeV. Chemical equi- 
librium requires that the chemical potential should satisfy 



M» -p p = p e = and p e = p^ 
Charge neutrality also requires 



(61) 



(62) 



that is, 



E 



7j 



2tt 2 A3 , 



EEV( £ i-«i) 2 -(l+27;4) 



2n 2 Xl 



P n=0 s? 
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Therefore, we can obtain the total mass-energy density of a 
npepir gas: 

P= ^2 Pj + m^n^. (64) 

j=n,p,e,/j. 

For a given p, we can determine all the quantities e e) e M , e p , e n , 
and n„ from Eqs. ( |6~ij ) - ( |64| ). Figure 7 shows the pion number 
density n n as a function of p for a given magnetic field. For 
l°g7c < 3, the magnetic field effect is negligible. 




15.1 15.3 15.5 15.7 15.9 16.1 



logp 



Fig. 7 The pion number density (in cm 4 ) as a function of p (in 
gem -3 ) for the log7 f = 4 (dashed line) and log7 f = 5 (dotted line). The 
solid line denotes the non-magnetic case (B=0). 

Similar to the muon appearance, we see that the threshold 
condition for 7r~ production is given as = 0. Finally, the pion 
production density in a magnetic field is given by 



47T 2 A 3 



(65) 



Figure 6 also shows the pion threshold density p* (B) as a func- 
tion of Y e . We also can see that p* (B) is not affected by mag- 
netic fields less than B ~ 10 17 G. However, above this field 
strength, the magnetic field effect on p* (B) is important. Note 
that since pions appear at the extremely high density of p ~ 10 15 
gem -3 , the AMM effect is small. 

From Eq. (|6l|), the pion momentum is zero (x^ = 0) because 
of the fact that in the ground state the condensed pions have 
zero kinetic energy. This forces the electron chemical potential 
e e to remain constant for p > p* . In consequence, the electron 
number density and pressure remain constant as p increases. 
Hence, increasing the pion density contributes to the total mass- 
energy density but not the pressure. As a result, for a given p, 
the pressure in the condensate phase is always less than in the 
non-condensate phase. Figure 8 shows the equation of state for 
an ideal magnetic npepir gas with pion condensation. We can 
see that magnetic fields reduce the pion condensation. How- 
ever, we still have a distinguishable pion condensate equation 
of state in strongly magnetized neutron stars. However, at this 
high density, the AMM does not contribute to the equation of 
state for the magnetized npepir gas although at j e > 10 5 the 
nucleon AMM changes somewhat the pion number density. For 



l°g7e < 3, the magnetic field effect is negligible and nearly the 
same as the non-magnetic npepir gas. Figure 9 shows the adia- 
batic index T as a function of p for a pion condensate equation 
of state. Here we also can see the oscillatory behavior of the 
adiabatic index. The inclusion of the AMM makes a secondary 
oscillation in the adiabatic index, but this secondary oscillation 
disappears gradually as density increases. 

7. DISCUSSION 

In this work, we have studied the nuclear equation of state 
for an ideal npe gas in a strong magnetic field. In particular, we 
have calculated the proton concentration, the threshold densi- 
ties for neutron, muon, and pion production and pion conden- 
sation in a strong magnetic field without and with the effect of 
the nucleon anomalous magnetic moments. In these calcula- 
tions, we have shown that the higher Landau levels are signif- 
icant at high density in spite of the existence of a very strong 
magnetic field. In particular, at high density, the proton concen- 
tration approaches the nonmagnetic limit. As a result, inverse (3 
decay is still suppressed in intense magnetic fields. Therefore, 
neutron-star rapid cooling is probably not affected by the direct 
URCA process which is enhanced in strong magnetic fields. In 
particular, we have obtained the neutron appearance density in 
a magnetic field when the nucleon AMM is included. We also 
show that the muon and pion appearance density are not af- 
fected by magnetic fields less than about B ~ 10 17 G. Finally, 
we here obtained an equation of state with a pion condensate 
in strong magnetic fields. Magnetic fields reduce the amount 
of pion condensation. However, we still have distinguishable 
effects of a pion condensate equation of state in strongly mag- 
netized neutron stars. In addition, we found the oscillatory be- 
havior of the adiabatic index in both strongly magnetized npe 
and npepiT gases at high density. Here we speculate that this 
behavior may also lead to an interior pulsational instability. 

It is generally accepted that neutrons and protons in a npe gas 
are superfluid. The c harged pion conden sate is also superfluid 
and superconductive ( |Migdal et al. 199C| ). This pion formation 
and condensation in dense nuclear matter wo uld have the sig- 
nificant consequence ( |Suh & Mathews 2000c ) that the equation 
of state would be softened. First of all, softening the equation of 
state reduces the maximum mass of the stars (B aym & Pethick 
1975). This softening effe ct with pion conden sation also leads 
to detectable predictions (Migdal et al. 1990). These are: (i) 
the enhanced rate of neutron star cooling via neutrinos, (ii) a 
possible phase transition of neutron stars to a superdense state; 
and (iii) sudden glitches in pulsar periods. In particular, if the 
pion condensation occurs in a strong magnetic field, it may sig- 
nificantly affect starquakes. 

According to the magnetar model by Duncan and Thomp- 
son (1992, 1995), SGRs are caused by starquakes in the outer 
solid crust of the magnetar. As a colossal magnetic field shifts, 
it strains the crust with huge magnetic forces and sometimes it 
cracks. When the crust snaps, it vibrates with seismic waves 
similar to those of an earthquake. However, in neutron stars 
they also produce a flash of soft gamma-rays. In addition, 
Cheng and Dai (1998) recently suggested that SGRs may be 
rapidly rotating magnetized strange stars with superconducting 
cores. 

Although such models can explai n some cruc ial features, 
there are still several unsettled issues (Liang 1995). Therefore, 
superconducting cores with a charged boson (pion and/or kaon) 
condensate in magnetars might be an alternative model to ex- 
plain the energy source of soft gamma-rays from magnetars. 




Fig. 9 Adiabatic index F as a function of p (in gem -3 ) for a pion con- 
densate equation of state for magnetic field strengths log7 e = 4 and 5. 
The thick and thin solid lines are for non-magnetic (B = 0) npe/j, and 
npefiir cases. 



FIG. 8 The equation of state for an ideal magnetic npefin gas with 
pion condensation. The thick and thin solid lines are non-magnetic 
cases (B = 0). The dashed and dotted lines are magnetic cases for 
log7 c = 4 and 5, respectively. 
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